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Abstract 

We study fluctuations of a Z7(l) gauge field on the worldvolume of Nf probe D7-branes in the background 
of N c black _D3-branes with a finite baryon density. The choice of mode corresponds to vector mesons in the 
dual gauge theory whose mass and lifetime can be determined by a study of the quasinormal modes. The 
speed of propagation of these holographic mesons at large momenta is studied from the dispersion relations 
of the quasinormal modes of the system. 



1 Introduction and results 



The AdS/CFT correspondence pQ provides a 
simple method to study strongly coupled gauge 
theories describing systems like the quark-gluon 
plasma formed at RHIC. Even in the absence of 
a concrete dual model describing all phenomena 
observed in QCD, holographic techniques can be 
used to study general features of systems with a 
large number of adjoint degrees of freedom. 

An example of such a setup is that of the back- 
ground geometry sourced by N c Z)3-branes (dual 
to J\f = 4 SYM theory) in which Nf < N c 
probe L>7-branes (describing quenched funda- 
mental J\f = 2 hypermultiplets) are embedded 
[2j. A thermalized plasma is modelled by the 
presence of a black-hole in the background AdS§ 
geometry. Finite baryon density, n_B, can be in- 
cluded in the model by turning on the temporal 
component of the abelian center of the U(Nf) 
gauge group present in the worldvolume of the 
probe branes. 

In the presence of baryon density the probe 
branes end on the black hole horizon [3J. This 
leads to finite lifetimes for the quarkonium 
bound states. The normal modes found in the 
absence of baryon density are unstabilized by its 



presence and can be described by a spectrum of 
quasinormal modes (QNMs) [4J. 

In [5], the low momentum dispersion rela- 
tion was studied by following the positions of 
the peaks in the spectral function caused by the 
QNMs. However, this method is not accurate 
for larger momentum due to numerical instabil- 
ities. Transforming the equation of motion for 
the QNMs into a Schrodinger form, the poten- 
tial, V(q), was studied. A similar study was car- 
ried out in [6]. A typical profile of this potential 
is found in figure [TJ 

In figure [T] an infinite barrier at the UV bound- 
ary (R* = 0) is observed. In the case of stable 
mesons (full line) there is also an infinite barrier 
at the tip of the embedded probe brane, thus giv- 
ing rise to a discrete spectrum (stable mesons). 
In [7] it was observed that these mesons prop- 
agate, at large momenta, at the local speed of 
light at the tip of the probe brane. 

When a finite baryon density is turned on, 
the IR barrier has a finite height (dashed line), 
and the range of the equation continues to the 
black hole horizon (i?* = — oo). At low mo- 
mentum, q, and frequency, u, one can observe 
quasiparticles in the spectral function, given by 
quasi-stable particles that tunnel the finite bar- 
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Figure 1: Typical profile of the Schrodinger potential 
for Minkowski (full line ) and black hole ( dashed line ) 
embeddings. The UV is at i?,, = and the horizon 
at i?* = — oo. In this case m — 4.53, q = 7.5 and 
ns = and 0.15 respectively. 



rier, leaking into the black hole. These are dual 
to relatively long-lived quasiparticles. If the fre- 
quency is raised far above the height of the bar- 
rier the particles lose their stability and become 
short-lived. 

When the momentum is increased the height 
of the barrier is lowered and a plateau at V oc q 2 
is formed. The exact momentum at which the 
barrier disappears defines a critical value, q cr it, 
which has been studied in detail in [5]. In this 
paper it was claimed that for q S> q cr i t the group 
velocity of the mesonic excitations should ap- 
proach the speed of light in the UV. In order 
to have such asymptotics the dispersion relation 
would need to change above q cr it- 

With an improvement in the numerical meth- 
ods we are able to study the large momentum 
regime, and find no signal of a change in be- 
haviour of the group velocity as suggested in [5] 
for q S> q C rit, see figure [2j 

Using a large-g approximation we show an- 
alytically, from the Schrodinger potential, that 
the group velocity for long-lived mesons (those 
with a large quark mass in the presence of a low 
baryon density) should be given approximately 
by the result in [71. i.e., the local speed of light 
of the same-m g Minkowski embedding (and thus 



for riB = 0) at the tip of the brane. 

2 Setup 

We consider the background geometry sourced 
by a stack of N c black D3-branes 

d* 2 = ^^p 2 (-£dt 2 + /d* 2 ) + ^d S | 6 , 

(i) 

where L = 4nrg s N c lg, T is the Hawking temper- 
ature, / = 1 - p- 4 , f = 1 + p- 4 and E 6 the 
6-dimensional euclidean space with isotropic ra- 
dius p. 

The Nf <C N c probe L>7-branes wrap an 
S 3 C S 5 of the background geometry. We de- 
scribe the Z)7-brane embedding by fixing one of 
the remaining angles = and expressing the 
other by a p-dependent function 6 = cos -1 x{p)- 

The action is given by 

Sdbi = N f T 7 f v /-detP[G + 27r^F], (2) 

with Tj the probe-brane tension, V the pull-back 
to the D7 worldvolume, G the 10-dimensional 
metric and F the field strength derived from the 
background gauge field A = A t {p) dt. 

The equation of motion for At can be solved 
in terms of a constant^ hb- Holographically, 
this constant corresponds to a baryon density, 
sourced by a chemical potential p = A t (pbou)- 
The equation of motion for x(p) is solved nu- 
merically, and its asymptotic expression in the 
UV gives the quark mass and condensate. 

The transverse vector mesons studied in this 
letter are given by the linearized fluctuation of 
the gauge field, A —* A + 5A(t, x 3 , p)dx 1 . See 
El E] for details on this construction. 

3 Analysis 

The pullback metric in the DBI action ([2D has 
the same holographic 4-dimensional manifold as 

1 7is coincides with the d of [5]- 
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(PQ), with the worldvolume coordinates given by 



ds D7 

L 2 



x 2 + p 2 x' 2 



L 2 



p\\ 



dp 2 +(l- X 2 )dfi| , 
(3) 

where dO| is the metric of a unit 3-sphere. 

Holography states that one can study the dual 
gauge theory by evaluating the on-shell action 
in the UV (p — > oo). An important point to 
note is that the metric ([3|) and the action ([2]) 
are invariant under p — > p^ 1 . This coordinate 
system covers the holographic space twice and 
does not describe the interior of the black-hole. 
This means that we can chose the range < p < 
1 for our studies. 

The choice of this domain greatly stabilizes 
the numerics as compared with the same calcu- 
lation in the 1 < p < oo domain. For large 
frequency and momentum, the standard checks 
of the numerical solution are satisfied. These 
include recovering the supersymmetric spectral 
function at large frequency [21 E], an exponen- 
tial decrease with frequency in the height of the 
peaks appearing in the spectral function [8], and 
the conservation of the Noether charge coming 
from the U(l) symmetry in the bilinear action 
for the fluctuation [TO] . 

To perform the numerical integration we chose 
the ingoing- wave boundary condition at the hori- 
zon. We then follow the poles of the retarded 
correlator by studying the values of the complex 
frequency u n , at which Dirichlet conditions in 
the UV {p — > 0) are satisfied [11] . The limit- 
ing group velocity is defined from this complex 
frequency as 



in) 



lim 

q— >oo 



dRe[u; n ] 
d q 



(4) 



Figure [2] shows the results for Re[u>i / (nT)] as 
a function of q/(irT) for cases I, II and III in 
table (3.19) of reference [5]. 

To understand why the limiting velocity does 
not approach the speed of light we express the 
equation of motion for the fluctuation in the 




Figure 2: Dispersion relation of the first mode for 
cases I (blue), II (pink) and III (yellow) in f^jj. 
The asymptotic velocities are respectively c g — 0.995, 
0.649 and 0.331. The lightcone lu = q serves to guide 
the eye. 



Schrodinger form [5] 

- d 2 R J(R*) +V(R*)ip(R* 



u 2 ^(R,) , (5) 



where i?* is a tortoise coordinate with the UV 
at i?* = and the horizon at i?* = — oo. The 
potential is split into two pieces V = Vq + V\q 2 , 
where < V\ < 1 (see [5] for the expressions 
of the functions). In the large momentum limit 
the piece proportional to q 2 is dominant every- 
where but close to the UV, where Vq diverges as 
(expressed in p coordinates) Vq ~ 

In the large momentum limit, for large m q and 
low baryon density ng, the potential takes the 
form given by the dashed line in figure [H with- 
out a barrier. The height of the potential can be 
approximated by V\(p e ff) q 2 , with p e jf an effec- 
tive radius giving the mean value of the height of 
the almost-flat plateau one can observe. Notice 
that Vi{p eff ) < 1. 

We can find an analytic result using a step 
function model with an infinite barrier at the 
origin. This was done for a different setup in 
[12] . from where we borrow the result 



V^pefrfqkkl + Oiq- 1 )] . (6) 



This gives precisely the dispersion relation found 
for Minkowski embeddings in [7]. In that case 
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the potential to take into consideration is the U- 
shaped one in figure [TJ which is approximated 
by a box of height V%(p e ff)q 2 . In the large q 
limit the solution coincides with ([6]) and p e ff is 
identified with the position of the tip of the probe 
brane. In [12] it was shown that an exponential 
tail on the left of the step or small bumps (such 
as those appearing for a momentum q < q cr it) do 
not affect this result, since these details enter at 
higher order q~ n in the large-g limit. 

When we consider vector mesons with short 
lifetimes (relatively low m q and/or large ng) the 
approximation taken for the Schrodinger poten- 
tial is not valid, and we expect the asymptotic 
group velocity to vary from the result in ([6]). In 
these cases the Schrodinger potential presents a 
ledge in the UV region, followed by a monoti- 
cally decreasing potential given by an approxi- 
mate constant gradient. At a given value of R* 
there is a kink and the value of the gradient 
changes dramatically, becoming approximately 
an exponential decay (see figure 20(b) in [5] for 
an example of the potential described above). 

Numerically it is observed that the value of 
Vi(R*) at the kink coincides with c 2 as obtained 
from the dispersion relation. This kink can be 
traced back to the qualitative change in be- 
haviour of the embedding profile x{p) ( see figure 
3 of [5]). With this result at hand, a study of 
the exact value p e ff(m q ,riB) at which the effec- 
tive velocity \fV\(p e ff) has to be evaluated to 
obtain the asymptotic group velocity is encour- 
aged. However, that calculation is beyond the 
scope of this letter. 

In this note we have given both numerical and 
analytical evidence to suggest that asymptotic 
velocities of quasinormal modes in the presence 
of baryon density do not necessarily approach 
the speed of light. 
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